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TECHNICAL MEMORANDUM X-64693 

TWO DIFFERENT APPROACHES FOR A CONTROL LAW 
OF SINGLE GIMBAL CONTROL MOMENT GYROS 

INTRODUCTION 

Double gimbal control moment gyros (DGCMG* s) are well suited for space 
applications as shown by W. B. Chubb and S . M. Seltzer (Reference 1) in their 
investigation of the Skylab Attitude and Pointing Control System. However , 
DGCMG' s have hardware problems since there are two gimbals and four bearings 
with each DGCMG. Thus, single gimbal control moment gyros (SGCMG* s) are 
now of increasing interest in the field of momentum exchange attitude control 
systems. The SGCMG' s are presently being considered for the High Energy 
Astrononomy Observatory (HEAO) and Large Space Telescope (LST) projects. 
Bendix, Sperry Rand, and TRW Systems have published reports relating to the 
use of SGCMG’s. (See References 2,3, and 4.) These reports, as well as simu- 
lations done by MSFC, show that problems occur from the singularities which 
appear in the SGCMG cluster. A singularity is present if the cluster cannot pro- 
duce a torque in a given direction in the three-dimensional space; then, at least in 
one direction, the spacecraft will not be controllable. The saturation singularity 
surface is given by the outer angular momentum envelope and can be avoided only 
by use of a desaturation system employing external torques. However, there are 
singularities within the angular momentum envelope which, theoretically, can be 
handled with desaturation by external torques. Thus, a specific control law for 
the SGCMG's is necessary to avoid the singularities. Some control laws have 
been investigated in the aforementioned reports with little success. In this report, 
two different approaches for the most general, nonlinear control law will be dis- 
cussed. The numerical computation of these control laws will not be given. A 
comparison between the two approaches will be made on a qualitative basis, since 
quantitative measures are a function of a particular CMG configuration. Thus, a 
comparison on a quantitative basis can only be made if a specified CMG cluster is 
given; however, some illustrations are included. 



MATHEMATICAL DESCRIPTION 


Let N be the number of single gimbal control moment gyros within the 

spacecraft. Let X = { , X 2 , X 3 } be an orthonormal frame fixed in the 

spacecraft and Z. = { Z* , Z„. , Z } an orthonormal frame fixed in the jth 
1 J 1 lj 2] 3j 

SGCMG, j = l( l)N. The frame Z^ is defined as follows (Fig. 1) : 




The angular momentum of the 
jth SGCMG can be described in tensor 
notation using the summation conven- 
tion 


Figure 1. SGCMG fixed 
frame. 


H.= V H..X.= H. Z. . , i= 1(1)3 , (1) 

J X ij l Zj ij ij v ’ 


where H. . and „.H.. are the coordinates of the angular momentum corre- 
X ij Zj ij 

sponding to frame X and Z^ , respectively. From the definition of frame Z_. 
for a completely balanced wheel, 

Zj H ij = I® h J ®] T - 

where h. is the jth SGCMG angular momentum. 


2 



It is practical to introduce a second spacecraft fixed frame Y = 

{ Yj , Y 2 , Yj which is correlated to the jth SGCMG in the following way: 

Yj. — axis is along the gimbal axis, 

Y ~ a? “ s is a ^ on S the zero gimbal angle direction, 

Y — axis completes the right-handed frame, 
dj 


Then the transformations between frame X , Y. , and Z. are charac- 


terized by 


3 


XY. 


XY. Y.Z. 


X. = J a., Y, . = 3 a.. 3 3 

i ik kj xk 


a M V k= 1(1)3 a = 1(1)3 • (3) 


XY. 

Here the matrix 3 a^ depends on the geometrical configuration of the jth 

XY. 

SGCMG within the cluster. Therefore, 3 a., is time-invariant. 

ik 


Y.Z. 

The matrix 3 3 a is an explicit function of the gimbal angle 6. of the 
!k j 


jth SGCMG: 


Y.Z. 
J 3 


1 

0 

0 

0 

eos6. 

- sin6. 


3 

3 

0 

sin5. 

cos6. 


a., = 

ik 


3 3 

Applying equation (4) to equation (2) one gets 
Y.Z. 

-,.II.. ~ 3 3 a. H.. =. [0 h.cosS. h.sinS.] 

Yj 13 ik Z. k] J J 3 3 


( 4 ) 


( 5 ) 


Here the theorem was used that vector coordinates are transformed just as the 
base vectors of the corresponding frames. 
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The angular momentum of the total CMG cluster will be evaluated in 
frame X which is the only frame used further, 


N XY. 


J 


X H i £ ~ik Yj"kj 

3=1 J 


a,,. ,,.H, . 


( 6 ) 


Using equation (5) , equation (6) can be rewritten using the convenient matrix 
notation, 


H = f(6) , H 

6 


'x 11 * ' x H * ’ x H ’’ 


K >1 » ^2 * • * • > ® 


(7) 


It results that the angular momentum of the CMG cluster is a nonlinear 
3x1 vector function f of the N x 1 gimbal angle vector 5 , From equation 
(7) the outer angular momentum envelope can be determined by variation of the 
gimbal angles. It has to be pointed out that the angular momentum envelope 
may have windows. Therefore, numerical investigations are recommended in 
each particular case. 


The torque that can be produced by the CMG cluster follows from the 
derivative of the angular momentum with respect to the inertial space. It 
yields 


T = H + w H , (8) 

where the 3 x 1 torque vector T and the 3x3 skew symmetric matrix a> 
appear. The matrix .£> corresponds to the angular velocity to of the frame X 
relative to the inertial frame. Usually the angular velocity is small and the 
second term in equation (8) is neglected. Thus, it remains 

# Qf • • 

T = H = _^6 = A6 , (9) 

where the 3 x N Jacobian matrix A of the nonlinear vector function appears. 
The CMG cluster is termed singular if the rank of matrix A is smaller than 
three; otherwise it may be termed regular: 


Rank A = 3 Regularity 
Rank A 3 Singularity 


( 10 ) 
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Usually singular surfaces can be found within the outer momentum 
envelope. If the angular momentum vector reaches a point on a singular surface, 
then the CMG cluster cannot produce a torque in an arbitrary direction. The 
resulting torque is constrained to a plane (Rank A = 2) or a line (Rank A = l) . 
This means that the controllability of the spacecraft is lost and arbitrary 
attitude errors may occur. Therefore, it is necessary to use at least three 
SGCMG' s, N S 3 , and it is very important to find a control law for the 
SGCMG’ s which guarantees the avoidance of singularities . 

For the formulation of a feasible control law some measure of the CMG 
cluster’ s regularity is needed. Such a measure can be found that is analogous 
to the degree of controllability introduced by P. C. Muller and H. I. Weber 
(Reference 5) using the quadratic symmetric matrix 

G - AA T = G T (11) 


One degree of regularity g is defined by the determinant of matrix G: 


g = det G = det AA^ 


( 12 ) 


Another degree of regularity is found from the eigenvalues A. of matrix G: 


g* = Min 


» j 1> 2, ... i N 


(13) 


Both degrees of regularity are positive for regular CMG clusters but they 
vanish in the case of singularities and depend on the gimbal angle vector 6. 


The numerical computation of equation ( 12) can be simplified using the 
Binet-Cauchy formula 


( N ) 

3 

i=l 


,N\ _ N(N-l)(N-2) 
u 3! 


N ^ 3 


(14) 
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where A . denotes the 


( N ) 

K 3 J 


distinct subdeterminants of order 3 extractable 


from the 3 x N matrix A. Thus, even in the case of numerous SGCMG' s 
there will be no essential difficulties in computation of the determinant 
equation ( 12) . 


In general a control law can be found by inversion of equation (7) or 
equation ( 9) . In the first approach, the gimbal angle vector 6 has to be 
controlled and the angular momentum vector H will be commanded. Using 
the second approach, the gimbal rate vector 6 is controlled and the torque 
vector T is commanded. 

GIMBAL ANGLE APPROACH FOR A CONTROL LAW 


The spacecraft will be controllable by the CMG cluster' s angular 
momentum if the Jacobian matrix A has full rank 3 for each angular momentum 
vector within the outer envelope. Optimal controllability will be given if the 
matrix A has an optimal regularity or, using equation ( 12) , if 


g( 6) = Max . (15) 

Thus, for each point within the angular momentum envelope, a gimbal angle 
vector 6 has to be found satisfying equation ( 15) . Since the gimbal angle 
vector is subject to equation (7) too, which can be interpreted as a constraint 
equation, necessary conditions for a local extremum can be found by the 
method of Lagrange multipliers. It yields (See Reference 6) 


H - f(<5) = 0 
T 

C - X A = 0 


(16) 


with the 1 x N Jacobian matrix C = 8 g/d 6, the 3 X N Jacobian matrix 
A = 9 f / 3 6 , the 3x1 vector X of the Lagrange multipliers, and N > 3. By 
equation (16) there are (3 + N) equations to solve for N gimbal angles and 
3 Lagrange multipliers. The solution of equation (16) will be a nonlinear 
N x 1 vector function, 
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6 = h(H) 


9 


(17) 


which can be named the inverse of equation (?) or the gimbal angle control 
law. 

The following cases now have to be discerned: 

1. h(H) and 3 h/3 H exist everywhere within the H-space of the 
envelope, 

2. h(H) exists and 3 h/3 H does not exist everywhere within the 
H-space of the envelope, 

3. h(H) does not exist everywhere within the H-space of the envelope. 

Case 1 . The commanded trajectory of the angular momentum vector 
can be achieved within the envelope by a smooth control of the gimbal angle 
vector. Specifically, by differentiation of equation ( 17) one obtains 

6 = || H = F(H)T , (18) 

where F is the N x 3 Jacobian matrix corresponding to equation ( 17) . 
Introducing equation (7) into equation ( 18) one obtains 

6 = F(6)T , (19) 

which defines an SGCMG control law without singularities within the 5 space 
of the cube: - ir — 6 — + tt, j = 1, 2, ... N > 3. 

Case 2 . The commanded angular momentum trajectory can be achieved 
within the envelope by a discontinuous control of the gimbal angle vector . This 
means that jumps occur in the gimbal angles at some surfaces within the 
envelope. During the time intervals needed for the jumps , the spacecraft will 
show some attitude errors. If the jump time interval approximates zero, the 
spacecraft errors vanish also,; however, infinite gimbal rates are necessary 
in the limiting ease , Such a jump phenomenon can be called internal desatu- 
ration because the CMG cluster is turned into a better regularity without any 
external torque. 
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Case 3 . The commanded trajectory of the angular momentum vector 
cannot always reach the outer envelope. There are holes within the 
envelope and, if the trajectory of the angular momentum vector is commanded 
to such a hole, external desaturation has to be provided by external torques. 

Programs for constraint parameter optimization problems can be used 
for the solution of equation ( 16) . Starting points may be the origin of the 
envelope and/or some known points on the envelope corresponding to the axis 
of the frame. 

GIMBAL RATE APPROACH FOR A CONTROL LAW 


The spacecraft can be controlled by the CMG cluster’ s torque if the 
Jacobian matrix A has the full rank 3 and especially if optimal regularity is 
required: 

g(5) = Max • (20) 

However, the torque equation, equation (9), with the gimbal rate vector 6 
as a controlled variable cannot be combined immediately with equation (20) 
for inversion. Thus, the general solution of linear equations first has to be 
treated, as by C. R. Rao andS. K. Mitra (Reference 7). 

A general solution of the consistent equation T = A 8 is 

6 = A~T + (E - A A)z , (21) 

where A is any N x 3 generalized inverse matrix of A , E is the N x N 
unit matrix, and z is an arbitrary N x 1 vector. 

Assuming the full rank 3 of matrix A the right inverse can be used 
as a generalized inverse matrix, 

A~ = a t (aa t ) -1 • ( 22 ) 
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Then, due to the properties of symmetric matrices, it yields 

E - A~A = E - AW)" 1 A = BB T , (23) 

T 

where B isan N x (N - 3) matrix. Further, Rank BB = N - 3 and AB s o. 
Since z is an arbitrary vector, it follows that 

T 

u = B z , (24) 


where u is an arbitrary (N - 3) x 1 vector, equally. From equation (21) 
to equation ( 24) , it follows that 


6 = aVV‘t + Bu 


(25) 


which is the inverse of equation (9) or the gimbal rate control law. The 
gimbal rate vector depends not only on the commanded torque T but also on 
an arbitrary vector u which produces no physical torque on the spacecraft. 
Thus, the vector u can be used for a redistribution of the gimbal angle vector. 
It is obvious from equation (.25) that the gimbal angles are not completely 
redistributable by u since T is an unknown, commanded function of time. 
However, if at least the degree of regularity is redistributable by u, 
then singularity avoidance seems to be possible. For the synthesis of a 
redistribution law, equation (20) has to be combined with equation (25) which 
can be done in two different ways: ( l) linear redistribution and (.2) nonlinear 
redistribution. 

Case 1 . From equation (20) it follows, examining the time history of 
the CMG cluster between zero angular momentum and outer angular momentum 
envelope, that 

4r s(6) = Max • (26) 


Since the degree of regularity g is a nonlinear function of the gimbal angle 
vector, it yields 
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g = 


y 


(27) 



C6 


with the lxN Jacobian matrix C known from equation (16). Using equation 

(25) , one obtains 

g - CA T (AA T ) _1 T + CBu . (28) 

Since the first term of equation (28) is commanded, only the second term of 
equation (28) can be used to redistribute the regularity and to satisfy equation 

(26) . There are two subcases possible: (l) D = CB = 0, locally unredistrib- 
utable and ( 2 ) D= CB + 0, redistributable. In case (l) the regularity g of 
the CMG cluster is locally unredistributable and therefore no redistribution law 
can be defined. In case (2) a suitable redistribution law is 


u = u q sgn d. ( 6) , (29) 


where d. are the elements of the .IX (N - 3) matrix D = CB depending on 

the gimbal angle vector 6. The redistribution gain u^ is arbitrary. However, 

by examining equation (28) it results that for an efficient redistribution, even 
in the neighborhood of a singularity g « 0, it is necessary that 



where u denotes the norm of the vector u . 

o o 


(30) 


Case 2 . For an optimal satisfaction of equation (20) , the absolute 
reachable regularity for each gimbal angle vector 6 within the gimbal angle 
cube has to be found. This can be done only in the limit case of zero torque 
T = 0 or infinite redistribution control ||u j| — » , because the commanded 
torque is an unknown function of time. Assuming zero torque T = 0, it 
remains from equation (25): 

6 = Bu . (31) 
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This nonlinear differential equation has to be solved for arbitrary redistribu- 
tion vectors u with all possible gimbal angle vectors 6 as initial conditions 
under simultaneous computation of the degree of regularity g, It is obvious 
that this can be accomplished, in general, only with considerable complexity. 
However, in the case of four SGCMG’ s, equation (3l) changes to two deter- 
mined equations 

6 = + Bu , 5 = - Bu , 

where u is a scalar. By solution of equation (32) there will be found, for 
each gimbal angle vector 6 within the admissible cube, the optimal sign of 
the redistribution control characterized by some number E. Thus, a redis- 
tribution law will be 


(32) 


u = sgn E(o) , u q -*• co . (33) 

However, in practice, infinite gimbal rates assumed for the synthesis of the 
redistribution law equation (33) are not feasible. Therefore, the optimal 
redistribution may be approximated only, since the additional influence of 
the commanded torque will change things continuously. 

The numerical computation of the linear control law equation (25) and 
the linear redistribution law equation ( 29) is simple. For the solution of the 
nonlinear differential equation (32) , resulting in the nonlinear redistribution 
law equation (.33) , presently available programs can be used. 


COMPARISON OF GIMBAL ANGLE AND 
GIMBAL RATE APPROACH 

A comparison of the gimbal angle and the gimbal rate approach is made 
in Table 1 assuming some simplifications. It is quite clear that the 
gimbal angle approach is much more suitable for the control law of the 
SGCMG’ s than the gimbal rate approach. The reason is that the fundamental 

T 

property of the CMG cluster, the regularity g - det AA , depends on the 
gimbal angles and not on the gimbal rates. Therefore, in the nonlinear 
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TABLE 1. COMPARISON OF GIMBAL ANGLE AND GIMBAL RATE APPROACH 
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gimbal rate approach, an integration is necessary to find a redistribution law. 
The linear gimbal rate approach, on the other hand, may be restricted due 
to the local information which results from the differentiation process 
characterizing the gimbal rate approach. In addition, a large torque command 
can overcome the provided redistribution and drive the CMG cluster to a small 
regularity. 

From the standpoint of technical realization, the linear gimbal rate 
approach is preferable because all computations can be done on line at the 
spacecraft. On the other side, the gimbal angle approach and the nonlinear 
gimbal rate approach need an intensive off-line computation and the resulting 
control laws have to be stored in the onboard computer in discrete form. 
However, the gimbal angle approach is defined over the three-dimensional 
angular momentum space, while the nonlinear gimbal rate approach is defined 
over the N-dimensional gimbal rate space. 


EXAMPLES 


For a regular four-CMG cluster, as shown in Figure 2, the mathemat- 
ical description will be given. For the illustration of the gimbal angle approach, 
an optimal gimbal angle vector will be determined for zero angular momentum 
and the linear gimbal rate approach will be treated in the neighborhood of a 
singularity. 


Let the four SGCMG’s have equal angular momentum in the correspond- 
ing frames Z.. Thus, from equation (2) one obtains 


Zj H.. = [0 h <)1 T , j = 1(1)4 


(34) 


Further, equation (5) results in 


_..tt . = h[0cos6 sin6 ] T 
Yj ij j j 


3 = 1(1)4 


(35) 
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( 36 ) 


( 37 ) 
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XY, 


a., = 
lk 


cos /3 

o 

sin/3 

*= ■ - 


0 

-1 

0 

> 

(38) 

sin/3 

0 

-cos/3 




and 


XY, 


a -i = 

lk 


cos/3 0 sin/3 

-sin/3 0 cos/3 

0-1 0 


(39) 


where /3 is a constant angle. From equation (6) and (7) one obtains 


H = 


sin/3 sinSj + sin/3 sinS 2 + sin/3 sinS 3 + . sin/3 sin5 4 

cosdj _ cos/3 sin S 2 - cos6 3 + cos/3 sinS 4 

cos/3 sindj + cos6 2 - cos/3 sin o 3 - cos5 4 


(40) 


and the Jacobian matrix, equation (9) , is 


A = 


sin/3 cosSj sin/3 cosS 2 sin/3 cos5 3 


-sinSj -cos/3 cos6 2 


sinOo 


cos/3 cosdj 


-sin6, -cos/3 cos6. 


sin/3 coso 4 
cos/3 cos6 4 


sin6 


4j 


(41) 


The minors A . of the matrix equation (41) are 


A. = sin/3 sinS. ^ sin(6. „ + 6. ) 

l 1+ 2 v i+1 i+3 7 

- sin/3 cos/3 cos 6. „ sin (5. .- 5 ) 

i+2 v i+1 i+3' 

+ 2sin/3 cos 2 /3 cos 6 cos 6 cos 6 

l+l i+2 i+3 

i = l( 1)4 and S 5 = 6 4 , 6 ( . = 6 2 , o 7 = 6 3 


(42) 
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For illustration of the gimbal angle approach, an optimal gimbal angle 
vector will be determined. In the case H = 0, there exists the particular 
solution 

61 = - 6 2 = §3 = - 64 = 0 • . (43) 

In regard to equation (43) one obtains from equation (42) 

= A = _ sin|3 sin0 sin 20 + 2 sinj3 cos 2 j3 cos 3 0 , 

i = 1(1)4 . (44) 


Thus, with equation (14) , the regularity of the CMG cluster for H = 0 is 


g = det AA T = Yj A . 2 = 4 A 2 (0) • (45) 

i=l 1 


The regularity, equation (45), is plotted for j3 = 45 deg in Figure 3. It is 
obvious that there are six singularities and two optimal gimbal angles. Thus, 
in this simple case one obtains 


6 = h(H= 0) = 0 . (46) 


Such an optimal gimbal angle vector has to be found for each angular momentum 
vector within the envelope. 

For illustration of the gimbal rate approach, a linearization in the 
neighborhood of the singularity <5j = S 3 - 0, S 2 = - 5 4 = - ir/2 is performed. 

At time t = t 0 let 
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6 


9 


(47) 


0 


®10 

7T 


• 

"2 + 9 ° 

• 

0 2O 


, 6 = 

• 

0 


0 3O 

7T 



_ 2 “ d °_ 


0 4O 


whcrti 0 O « 1 and 0 .^ , j - l( 1)4 , is arbitrary. Then, one obtains from 

equation (41) the following Jacobian matrix neglecting all terms of third and 
higher order: 


sin/3 

sin/30 o 

sin/3 

sin/30, 

0 

-cos|30 o 

0 

cos/30, 

cos/3 

1+ | 0 O 2 

-cos/3 



with the regularity 

g = det AA^ = 8 sin 2 /3cos 2 /3 (l ■+ eos 2 /3)0 o 2 


(48) 


(49) 


Further, it yields 

sin/3 cos 2 /3(l + cos 2 /3)0 o 2 0 sin 2 /3cos 3 i30 o 2 

0 -sin 2 /3 cos/3 ( 1 + cos 2 /3)0 o 2 sin 2 /3cos 2 /30 o 2 

sin/3cos 2 j3(l + cos 2 /3)0 o 2 0 -sin 2 /3cos 3 /30 o 2 

0 sin 2 /3cos/3(l + cos 2 /3)0 o sin 2 /3cos 2 /30 o 2 

( 50 ) 


a t (aa t ) 


-i 


4 

g 
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and 

-sin/3 cos/3 0 O - sin/3 cos 2 /30 o 2 
sin/3 cos z /3 0 O 

sin/3 cos/3 0 O - sin/3 cos 2 /30 o 2 
sin/3 cos 2 /3 0 o 

The matrix C according to equation (27) is 

C = 8 sin 2 /3cos 2 /3(l + cos 2 /3 )0 O [0, 1, 0, - 1] 



( 51 ) 


(52) 


and the matrix D, in this example a scalar, is 

D = CB = 0 . (53) 

Thus, it results that the CMG cluster is locally unredistributable even though 
the regularity, equation ( 49) , does not vanish. 

CONCLUSION 


Two different approaches for a control law of single gimbal control 
moment gyros have been presented but the gimbal angle approach seems to be 
better suited for a control law than the gimbal rate approach because of the 
singularity problem. The gimbal angle approach requires off-line compu- 
tation of suitable gimbal angles for each point within the three-dimensional 
angular momentum envelope and the commanded variable is the angular 
momentum itself. The linear gimbal rate approach can be implemented 
by online computation and the commanded variable will be the torque. The 
nonlinear gimbal rate approach also requires off-line computation of suitable 
redistribution controls for each point within the N-dimensional gimbal angle 
cube. However, a general statement on the controllability of spacecrafts 

with an SGCMG cluster has not been found. Only detailed numerical compu- 
tations will answer this question. 
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